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Abstract 


Flood is one of the natural disasters that can threaten economic, social and environmental 
sustainability activities. Therefore, identifying flood-prone areas in the basins is essential. The 
aim of this research was to prioritize and analyze flooding potential in Ghamsar watershed using 
multi-criteria decision making models and fuzzy logic. For this purpose, first, 14 effective 
criteria on basin flooding (including precipitation, elevation, land use, soil texture, lithology, 
distance to river, slope, drainage density, stream order, Topographic Wetness Index, Stream 
Power Index, flow accumulation, plan curvature, profile curvature) were identified. Rhen, 
pairwise comparisons were made between the criteria using the AHP method in Expert choice 
11 software and the final weight of each criterion was obtained. Moreover, the criteria were 
fuzzified using linear and large fuzzy functions, and finally, the final weights obtained for each 
criterion were applied in the respective layers and the final flood risk map of the basin was 
prepared. The results of prioritizing quantitative and qualitative criteria based on the opinion of 
watershed management experts using the AHP method showed that among the 14 proposed 
criteria, the precipitation criterion with a final weight of 0.229 has the most impact on the flood 
risk of the basin. Moreover, the criteria of heavy to very heavy soil texture (0.499), Plvav stone 
unit (0.252), water areas and bare lands with final weights of 0.345 and 0.225, respectively, 
have a stronger role in the flood risk of the basin. The results of the final flood risk map of the 
basin showed that about 30.841, 27.056 and 12.406 percent of the total area of the basin are in 
the medium, high and very high flood risk categories, respectively, and these areas are located 
mostly in the central part and along the southeast of the basin. Therefore, knowing flood 
potential of the basin can be effective in formulating crisis management plans when faced with 
floods. 
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جغرافیا و مخاطرات محیطی. سال دوازدهم. شمارة چهل و هشتم. زمستان ۰۱۶۰۲ صص ۱۳۹-۱۵۹ 


اولویت‌بندی و تحلیل فضایی پتانسیل سیل‌خیزی مبتنی بر رویکرد FUZZY-AHP‏ 
(مطالعه موردی: حوزه آبخیز قمصر) 
دانیال صیاد- دانش آموخته دکتری علوم و مهندسی آبخیزداری- آب» دانشکده منابع طبیعی و علوم زمین, دانشگاه کاشان 
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تاریخ دریافت: ۱٤۰۰/۱۲/۱٩‏ تاریخ بازنگری: ۱۶۰۱/۷/۲ تاریخ تصویب: ۱۶۰۱/۳/۱۷ 


حکده 


سیل یکی از بلایای طبیعی اسست که می‌تواند فعالیت‌های اقتصادی. اجتماعی و پایداری محیط را تهدید 
کند؛ بنابراین شناسایی مناطق مستعد سیل‌خیز حوزه آبخیز امری ضروری در این زمینه است. هدف پژوهش 
حاضر اولویت‌بندی و تحلیل فضایی سیل‌خیزی حوزه آبخیز قمصر با استفاده از مدل‌های تصمیم‌گیری 
چند معیاره و منطق فازی است که با رویکرد توصیفی — تحلیلی و با روش FUZZY-AHP‏ انجام‌شده 
است. برای این منظور ابتدا ۱۶ معیار موثر بر سیل‌خیزی حوضه (شامل باران» ارتفاع کاربری اراضی» بافت 
خاک» سنگ‌شناسی, فاصله از رودخانه» شیب. تراکم زهکشی, رتبه آبراهه» شاخص رطوبت توپوگرافی 
شاخص توان آبراهه» تجمع جریان, انحنای سطح» انحنای طولی) شناسایی شد و سپس در نرم‌افزار 
Expert choice 1‏ بین معیارها با روش AHP‏ مقایسه‌های زوجی انجام گردید و وزن نهایی هر معیار به 
دشت آمك. در ادامه» معیارها با استفاده از توابع فازی large linear‏ فازی‌سازی شدند و درنهایت 
وزن‌های نهایی به‌دست‌آمده برای هر معیار در لایه‌های مربوطه اعمال شد و نقشه نهایی سیل‌خیزی حوضه 
تهیه گردید. نتایج حاصل از اولویت‌بندی معیارهای کمی و کیفی بر اساس نظر کارشناسان آبخیزداری با 
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۱۳۰ جغرافیا و مخاطرات محبطی شمارة چهارم 


روش AHP‏ نشان داد از بین ۱۶ معیار پیشسنهادی, معیار باران با وزن نهایی ۰/۲۲۹ بیش‌ترین تأثیر را در 
سیل‌خیزی حوضه دارد. همچنین معیارهای بافت خاک سنگین تا خیلی سنگین (۰/1۹۹) واحد سنگی 
Plvav‏ ۰۰/۲۵۲ کاربری‌های پهنه‌های آبی و اراضی لخت به ترتیب با وزن نهایی ۰/۳۶۵ و ۰/۲۲۵ QAI‏ 
پررنگ‌تری در سیل‌خیزی حوضه دارند. نتایج حاصل از نقشه نهایی سیل‌خیزی حوضه نشان داد که حدود 
۱ ۷/۰۵ و ۱۳/۰۲ درصد از مساحت کل حوضه به ترتیب در طبقه سیل‌خیزی متوسط زياد و 
خیلی زیاد قرار دارند که این مناطق بیشتر در بخش مرکزی و در امتداد جنوب شرقی حوضه واقع شده‌اند؛ 
بنابراین آگاهی از پتانسیل سیل‌خیزی حوضه می‌تواند در جهت تدوین برنامه‌های مدیریت بحران در هنگام 
مواجه‌شدن با سیلاب موثر واقع شود. 

کلید واژه‌ها: توابع فازی, مناطق مستعد سیل خیز مناطق نیمەخشک. نرم‌افزار Expert choice‏ 


۱- مقدمه 

سیل را می‌توان جریان‌هایی که فراتر از ظرفیت حمل‌ونقل کانال و دریاچه اسست. تعریف نمود که می‌تولند به 
دلایل طبیعی و مصنوعی مانند بارندگی شدید. شکست ناگهانی سد در حوضه‌ها و غیره باشد (تنزین و باسکار" 
۰ سیل با آسیب‌های وارده به ساختمان‌هاء اقتصاد و فعالیت‌های اجتماعی پایداری محیط را تهدید می‌کند که 
منجر به از دست رفتن زندگی می‌شود (اکمکچی اوغلو" و همکاران» ۲۰۲۱). ازاین‌جهت سیل به فاجعه‌ای تبدیل‌شده 
است که اثرات زیان‌بار آن از سایر بلایای طبیعی ازجمله زمین‌لرزه بیش‌تر اسست (کای " و همکاران, ۲۰۱۹). طبق 
بررسی اولیه پایگاه داده فاجعه بین‌المللی از سال ۱۹۰۰ تا ۲۰۱۹ حدود هفت میلیون نفر در اثر سیل کشته شده‌اند و 
همچنین بیش از ۷۰۰ میلیارد دلار خسارت در سراسر جهان برجا گذاشته است Sly)‏ و همکاران ۲۰۲۰). با توجه 
به مطالب بیان شده بیش‌تر آسیب‌های ناشی از بلایای سیل در مقیاس محلی رخ می‌دهد که بر زندگی مردم؛ 
زیرساخت‌ها؛ کشاورزی و ساختار اقتصادی اثر می‌گذارد که درنهایت» توسعه اقتصادی محلی را محدود می‌کند و 
ازاین جهت. ارزیابی خطر سیل برای کاهش خسارات ناشی از سیل امری ضروری است (گوا* و همکاران» QE‏ 
امروزه چندین روش برای تخمین حساسیت سیل در منطقه وجود دارد که هر رویکرد. جنبه‌های منفی و مثبت خود 
را دارد. در این راستا می‌توان به رویکردهای تصمیم‌گیری چند معیاره بر پایه تحلیل سلسله مراتبی (AHP)‏ منطق 
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فازی CFL)‏ رگرسیون لجستیک LR)‏ درخت تصمیم MDT)‏ یادگیری ماشین EML)‏ شبکه عصبی مصنوعی 
(ANN)‏ اشاره کرد که به‌طور گسترده موردپذیرش محققین قرارگرفته است (مالیک" و همکاران» ۲۰۲۰). 
تصمیم‌گیری چند معیاره" بر پلیه مجموعه‌ای از روش‌ها و فرآیندهایی است که برای ارزیابی و انتخحاب بهترین 
گزینه‌های تصمیم بر اساس معیارهای متعدد استوار است (فرناندز" و همکاران, ۲۰۱7). فرآیند سلسله مراتبی یک 
تحلیل ساختاری است که تعصب در تصمیم‌گیری را کاهش می‌دهد و می‌تواند ابزار مناسبی در برابر تصمیمات 
پیچیده باشد (کولیموشی* و همکاران. ۲۰۲۱). همچنین منطق فازی یکی از رویکردهای مرسوم برای تعیین حساسیت 
سیل در منطقه است که برای مقابله با مشکلات پیچیده کارآمد است (مالیک و همکاران. ۲۰۲۰). امروزه مفهوم منطق 
فازی در علوم محیطی به یک روش آماری پیشرفته و مهم تبدیل شده است (سینشاو " و همکاران, ۲۰۲۱). 
رویکرد FUZZY-AHP‏ یک روش بهبودیافته برای تنظیم وزن‌های حاصل از AHP‏ بر اساس مجموعه فازی و منطق 
فازی است که خطر ابتلا به سیل را نسبت به روش AHP‏ معقول OUS‏ می‌دهد (وانگ و همکاران. ۲۰۲۱). ازاین جهت 
در پژوهش حاضر از FUZZY-AHP» Sus)‏ برای اولویت‌بندی تحلیل فضایی سیل‌خیزی حوزه آبخیز واقع در 
منطقه نیمهعشک (حوضه قمصر کاشان) استفاده شد. تاکنون مطالعات متعددی در زمینه شناسایی مناطق مستعد 
سیل‌خیزی با تحلیل سلسله مراتبی و منطق فازی انجام شده است که ازجمله این پژوهش‌ها می‌توان به تحقیق‌های 
زیر اشاره کرد: مصابی و مکونی (THT)‏ طی تحقیقی در تانزانیا برای agi‏ نقشه حساسیت سیل از تجزیه‌وتحلیل 
تصمیم‌گیری چند معیاره و سامانه اطلاعات جغرافیایی استفاده کردند. نتیجه مطالعات آن‌ها نشان داد که به ترتیب 
حدود ٤٤/٩۷‏ و ۵۹/۷۳ درصد در محدود احتمال خطر سیل زیاد تا خیلی زیاد و خیلی کم تا متوسط قرار دارند. 
Tbe‏ و همکاران (۲۰۱۹) در پژوهشی باهدف منطقه‌بندی خطر سیل با استفاده از تحلیل سلسله مراتبی و تثوری 
فازی بیان داشتند. احتمال خطر از شمال غربی به جنوب شرقی به‌تدریح افزایش‌یافته است. همچنین این محققین بیان 
داشتند حساس‌ترین مناطقء مناطقی با توپوگرافی نسبتاً کم بودند. یودیینگ"' و همکاران (۲۰۱۹) به بررسی خطر 
سیل با استفاده از فرآیند تحلیل سلسله مراتبی و تحلیل فازی در تایلند پرداختند. نتیجه پژوهش آن‌ها OUS‏ داد بیش‌تر 
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مناطق نزدیک به شبکه زهکشی در be tol‏ خیلی زیاد قرار دارند. زیاری و همکاران (۱2۰۰) به‌منظور پهنه‌بندی 
پتانسیل سیل‌خیزی با روش تحلیل سلسله مراتبی و منطق فازی در شهر ايلام به این ame‏ دست يافتند که نواحی 
شمال و شمال غربی بیش‌ترین ظرفیت در برابر خطر سیل را دارند. محمدپور و وفایی‌نژاد (۱۳۹۹)» طی پژوهشی 
به‌منظور اولویت‌بندی و شناسایی عوامل موثر بر رخداد سیل در حوضه تجن با استفاده از تحلیل سلسله مراتبی و 
منطق فازی به این نتیجه رسیدند که مهم‌ترین تمرکز مناطق حساس سيل در محدوده رودخانه اصلی و محدوده دلتای 
رودخانه در مصب است. نجفی و کریمی کرد آبادی YA)‏ در پژوهش خود با استفاده از مدل ترکیبی AHP-‏ 
FUZZY‏ به ارزیابی و پهنه‌بندی خطر سیلاب در منطقه یک کلان‌شهر تهران پرداختند و به این amem‏ دست پافتند 
پهنه‌های با خطر خیلی کم. خطر کم و متوسط در شمال و مرکز محدوده موردمطالعه قرار دارند. همچنین این محققین 
بیان داشتند پهنه‌های با خطر خیلی زیاد و زیاد در محدوده شهری منطقه و منطبق بر نقاط خروجی حوضه‌ها است. 
صالحی و همکاران (۱۳۹۲) با استفاده از فرآیند تحلیل سلسله مراتبی فازی به پهنه‌بندی خطر سیلاب شهری در 
تهران پرداختند و به این نتیجه دست یافتند نواحی مرکزی تا جنوب شهر و همچنین آبراهه‌ها و حرایم آن‌ها بیش‌ترین 
پتانسیل ازلحاظ سیل‌خیزی را دارند. علی‌رغم مطالعات فراوان در زمینه پهنه‌بندی خطر سیلاب در مناطق مختلف. در 
کم‌تر پژوهشی به بحث سیل‌خیزی در مناطق نیمه‌خشک پرداخته شده است؛ بنابراین پژوهش حاضر با هدف 
اولویت‌بندی و ناحیه‌بندی مناطق مستعد سیل‌خیز با انتخاب معیارهای گسترده و موثر بر سیل‌خیزی حوضه ازجمله 
باران. سنگ‌شناسی» کاربری اراضی (تصاویر ستتینل ۲۲ ۰ شیب ارتفاع بافت خاک تراکم زهکشی فاصله تا 
آبراهه رتبه آبراهه‌هاء تجمع جریان. شاخص رطوبت توپوگرافی . شاخص توان آبراهه " انحنای سطح و انحنای 
طولی" در تعیین پتانسیل سیل‌خیزی حوزه آبخیز واقع در مناطق نیمه‌حشک با استفاده از رویکرد تصمیم‌گیری چند 
معیاره (MCDA)‏ و منطق فازی (FUZZY LOGIC)‏ انجام شده است. 


۲- مواد و روش 
۱-۷- منطقه مو ردمطالمه 


حوزه آبخیز قمصر با مساحت حدود۱ ۲۶ کیلومترمربع در محدوده بین طول جغرافیایی YA yy" boy? ۱۸ ol‏ 
019 شرفی و عرض جغرافیایی YY? YA! A"‏ تا ov YA"‏ ۳ شمالی قرار دارد J=)‏ ۱). بر اساس آمار ۲ ایستگاه 
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2 Topographic Wetness Index 
3 Stream Power Index 

4 Plan Curvature 
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شکل ۱- موقعیت منطقه موردمطالعه 
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پژوهش حاضر از نوع توصیفی- تحلیلی و با روش FUZZY-AHP‏ و در نرم‌افزارهای 11 Expert choice‏ و 
Arc Gis10.5‏ انجام شده است. برای این منظور. ابتدا داده‌های موردنیاز ازجمله نقشه طبقات ارتفاعی منطقه (DEM)‏ 
با اندازه سلولی ۱۲/۵ × ۱۲/۵ agi‏ شد. نقشه کاربری اراضی با اندازه سلولی ۱۰ × pele‏ برای سال ۲۰۲۰ از تصاویر 


ستتینل wo ang Y‏ در جدول(۱)» منابع op‏ داده‌های ad ol‏ آورده Aa ados‏ همچنین نقشه بافت خاک منطقه از 


زوجی وارد شدند و با قضاوت کارشناسان (اساتید گروه آبخیزداری و دانشجویان دکتری آبخیزداری) » وزن ale‏ هر 
لایه و زیرمعیارهای لایه‌های کیفی (نقشه بافت خاک کاربری اراضی و واحد سنگ‌شناسی) با در نظر گرفتن نرخ 


سازگاری در قضاوت‌ها Jus Oe‏ در پژوهش حاضر برای مقایسه‌های زوجی از مقیاس عددی ۱ 5 q‏ استفاده شد 


۱۳۴ جغرافیا و مخاطرات محبطی شمارة چهارم 


که مقایسه Y‏ بیانگر ترجیح برابر و مقایسه ٩‏ بیانگر ترجیح مطلق در مقایسه‌های زوجی است W)‏ و بالوگون؛ ۲۰۲۰). 
همچنین برای ارزیابی قضاوت‌های انجام‌شده (مقایسه‌های زوجی) از شاخص ناسازگاری استفاده شد. سپس 
معیارهای مؤثر بر تحلیل فضایی سیل‌خیزی حوضه قمصر اولویت‌بندی شدند. درنهایت لایه‌های انتخابی (معیارها) با 
استفاده از ابزار Fuzzy Membership‏ در نرم‌افزار 10.5 Arc Gis‏ استانداردسازی (فازی) شده و در درجه عضویت 
صفر تا یک قرار گرفتند. ازآنجایی که مقادیر اندازه‌گیری شده از معیارها ممکن است در واحد و دامنه‌های متفاوتی به 
ثبت رسیده باشند ALL‏ تمام bay‏ در دامنه همسانی ۰ تا ۱ جا به جا شوند تا استانداردسازی شوند. به‌پیان‌دیگر تابع 
عضویت ply‏ ۱ بیانگر عضویت کامل در مجموعه و تابع عضویت صفر Sly‏ عضویت نداشتن کامل در مجموعه 
فازی است (زیاری و همکاران» ۱۶۰۰). همچنین فازی سازی نقشه‌های aS‏ استفاده از وزن نهایی طبقات آن‌ها در 
نرم‌افزار 11 Expert choice‏ با روش تحلیل (AHP) vu alala‏ و قضاوت کارشناسی تعیین و با استفاده از تابع 
فازی large‏ فازی سازی شدند. در ادامه برای تحلیل فضایی سیل‌خیزی با GA‏ وزن‌های به‌دست‌آمده از روش تحلیل 
سلسله مراتبی (AHP)‏ در نقشه‌های فازی شده هر (Les) LY‏ نقشه سیل ‌گیری حوضه 4 دست EX‏ درنهایت 
به‌منظور تحلیل فضایی‌سیل‌خیزی حوضه نتایج در پنج گروه خیلی کم. کم متوسط. زیاد و خیلی زياد طبقه‌بندی 
شدند yl»)‏ ۲ و همکاران Q4‏ در Ys‏ مراحل انجام پژوهش نشان داده شده است. 


حدول ۱- منابع 4d sl‏ تهیه داده‌ها 


ردیف داده دقت مکانی منبع 
۱ مدل رقومی ارتفاعی ۵ x‏ ۱/۵ متر https://search.asf.alaska.edu‏ 
۲ نقشه کاربری اراضی ۰ ۱ متر Esri 2020 Land Cover‏ 
۳ نقشه خاک شهرداری کاشان 
£ داده‌های باران‌سنجی شرکت آب منطقه‌ای اصفهان 


1 Tella & Balogun 
2 Dano 
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م تحلیل فضایی سیلخیزی حوضه 


v 


ے معیارهای موثر بر سیلخیزی حوضه 


v 


} 
ابه های کیفی ^4 فازی سازی لابه ها “ye‏ های کمی 


شکل ۲- مراحل انجام پژوهش 


۱-۲-۲-تاثیر هر معیار بر سیل‌خیزی حوزه آبخیز و فازی سازی لایه‌ها 

در پژوهش حاضر انتدا با بررسی در پژوهش‌های مختلف (تلا و evo «o XS JL.‏ داس T1۸ e‏ مالیک و 
همکاران» DERE‏ مصابی و مکونی. eY Y)‏ سویین" و (Yeye «UA‏ معبارهای مور بر سیلاب شناسایی dico‏ 

بارندگی 

SL‏ مهم‌ترین عامل در ایجاد سیل است (بوامرانه " و همکاران ۲۰۲۰). هرچه شدت و فراوانی بارندگی 
بیش‌تر باشد احتمال وقوع سیلاب بیشتر می‌شود (تلا و بالوگون. ۲۰۲۰). در این تحقیق, برای فازی‌سازی نقشه پهنه 
باران از تابع حطی * - افزاینده استفاده شد. همچنین نقشه پهنه باران با استفاده از میانگین بارندگی سالائه ۳۷ ساله) 
سه ایستگاه بن رود فین و برزک (جدول ۲) به روش IDW‏ (معکوس فاصله)* در ag Arc Gisl0.5‏ شد. 


1 Das 

2 Swain 

3 Bouamrane 

4 linear 

5 Inverse Distance Weighting 


شمارة چهارم 


جغرافیا و مخاطرات محیطی 


lY? 


جدول ۲- موقعیت ایستگاه‌های باران‌سنجی 


ایستگاه طول جغرافیایی عرض جغرافیایی ارتفاع 
بن رود ۳ ۷ yr? ££ £v" o‏ ۱۹۰ 
برزک YY? £V. YA" oy? ۱۳ yg"‏ ۸۰ 
فين YY? ov va" oy? vv vv‏ ۱۰ 


بر سیل‌خیزی حوضه. در این پژوهش برای فازی‌سازی این نقشه از تابع ئ استفاده شد. در حدول ۲ مشخصات 


سنگ‌شناسی حوضه ارائه شده است. 


حدول ۳- مشخصات سنگ‌شناسی منطقه موردمطالعه 


واحد سنگی مشخصه m‏ ساحت (km?)‏ 
Etliv‏ توف آتش‌فشانی آندزیتی ائوسن AVY‏ 
Klaml‏ سنگ‌آهک خاکستری با بستر ضخیم تا Sal Se‏ اوربیتولینا کرتاسه ۱/۳ 
Oüv‏ دیوریت الیگوسن 1۷ 
OMI‏ سنگ‌آهک مرجانی الیگوسن- میوسن | ۵۰/۰۳ 
OMsm‏ | سنگآهک» مارن مارن گچی, ماسه شنی و ماسه‌سنگ (سازند قم) | الیگوسن میوسن Ví‏ 
OMviv‏ آتش‌فشان آندزیتی الیگوشن vv‏ 
PLvav‏ آتش‌فشان‌های ریولیتی تا ریوداسیتی پلیوسن Utv‏ 
Q2‏ نهشته‌های تراسی کواترنر NY‏ 
 URig‏ | مارن قرمن مارن گچی» ماسه‌سنگ و کنگلومرا (سازند قرمز بالایی) میوسن 0 


] Patrikaki 
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کاربری اراضی 

کاربری اراضی عامل مهمی در مطالعات خطر سیل است و می‌تواند بر اجزای فرآیندهای هیدرولوژیکی حوضه 
مانند رواناب نفوذ و تبخیر و تعرق تأثیر بگذارد (سوئیسی! و همکاران ۲۰۲۰). ازآنجایی که نقشه کاربری اراضی 
ماهیت کف aly ute‏ فازی‌سازی abd‏ کاربری اراضی کی زرهق ساقس اعدا با اسفاده از فقضاوت کارشناسی و 
مقایسه‌های زوجی بین هر نوع کاربری با فرآیند تحلیل سلسله مراتبی. وزن نهایی هر نوع کاربری تعیین شد و سپس 
این اوزان در نقشه کاربری اراضی اعمال شد و درنهایت به نقشه رستری تبدیل گردید. سپس کاربری‌ها با توجه به 
طبقه ستتینل ۲ در ۸ طبقه به‌صورت جدول ٤‏ طبقه‌بندی شدند. برای فازی‌سازی نقشه کاربری اراضی نیز از 
تابع Large‏ در Arc Gis10.5‏ استفاده شد. 


جدول -٤‏ طبقه‌بندی کاربری اراضی بر اساس تصاویر سنتینل ۲ 


(km?) کاربری اراضی مساحت‎ 
4 Water 
JAM Trees 
SESS Grass 
st Flooded vegetation 
۳۵ Crops 
3/110 Scrub/shrub 
VIVEAN Built Area 
YE/VOA Bare ground 


: t 


با کاهش شیب (سطوح مسطح) سرعت جریان آب کاهش يافته و باعث تجمع آب می‌شود. ازاین جهت مناطق 
مسطح نسبت به مناطق شیبدار از حطر ابتلا به سیل‌خیزی بالاتری برخوردار هستند (ریمبا" و همکاران. ۲۰۱۷). در 
پژوهش حاضر نقشه شیب از روی نقشه مدل رقومی ارتفاع منطقه " در 0910.5 Arc‏ تهیه گردید. همچنین برای 
فازی‌سازی نقشه شیب از تابع خطی — کاهنده استفاده شد. 


1 Souissi 
2 Rimba 
3 Digital Elevation Model 


۱۳۸ جغرافیا و مخاطرات محیطی oles‏ چهارم 


خاک 

بافت خاک با توجه نقش آن در نفوذ آب به داخل خاک بر وقوع سیل در منطقه موثر است؛ به‌طوری‌که خاک‌های 
درشت‌دانه (بافت سبک)» سیل‌گیری کم‌تری نسبت به خاک‌های ریزدانه (بافت سبک) دارند (ابراهیم بتیس و احمد؛ 
۲ ازآنجایی که لایه بافت خاک یک لایه AS‏ محسوب می‌شود برای فازی‌سازی لایه خاک از تابع Large‏ در 
Arc 55‏ استفاده شد. 

ارتفاع 

آب به دلیل نیروی گرانش به‌سرعت از مناطق بالادست و ارتفاعات بالاتر به سمت مناطق پایین‌دست حرکت 
می‌کند و با رسیدن آب به مناطق پایین‌دست (مناطق کم ارتفاع) گسترش می‌یابد و باعث سیل‌خیزی می‌شود (داس؛ 
۰ در این پژوهش با توجه به تأثیر عامل ارتفاع در سیل‌خیزی حوضه برای فازی‌سازی AY‏ ارتفاع از تابع Be‏ 
- کاهنده استفاده شد. 

تراکم زهکشی 

یکی از عوامل اصلی که باعث وقوع سیل می‌شود تراکم زهکشی است که با بالا رفتن تراکم زهکشی میزان خطر 
ناشی از سیلاب افزایش می‌یابد (سوئیسی و همکاران» ۲۰۲۰). نقشه تراکم زهکشی در این پژوهش در 0510.5 Arc‏ 
و با ابزار Line distance‏ تولید شد. سپس برای فازی‌سازی نقشه مذکور از تابع خطی - افزاینده استفاده شد. 

فاصله از آپراهه 

مناطق نزدیک به شبکه زهکشی رودخانه به‌شدت تحت تأثیر خطر سیل قرار دارند. درحالی‌که تأثیر این پارامتر بر 
سیل‌خیزی حوضه با فاصله گرفتن از شبکه زهکشی کاهش می‌یابد (فرناندز و لوتز ۲۰۱۰۵). ازاین‌جهت مناطق 
نزدیک به رودخانه پتانسیل سیل‌خیزی بالایی دارند و برای فازی سازی نقشه فاصله از آبراهه از تابع Jas‏ — کاهنده 
استفاده شد. همچنین نقشه فاصله از آبراهه با استفاده از ابزار Euclidean Distance‏ در Arc Gis10.5‏ تولید شد. 

رتبه‌بندی آبراهه 

با بزرگ‌تر شدن رتبه آبراهه شدت» سرعت و دبی جریان افزايش می‌یابد (کاستا؛ ۱۹۸۷ بنابراین آبراهه‌های با 
رتبه بالاتر از پتانسیل ایجاد سیل بالاتری برخوردارند. در این مطالعه نقشه رتبه آبراهه با استفاده از ابزار Hydrology‏ و 
با روش استرالر در نرم‌افزار Arc Gis10.5‏ تهیه گردید و برای فازی‌سازی این og «M‏ برای فازی‌سازی «Y‏ رتبه 
aal‏ از تابع خطی - افزاینده استفاده شد. 


1 Ibrahim-Bathis & Ahmad 
2 Fernandez & Lutz 
3 Costa 
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تجمع جریان 
تجمع جریان معیار مهمی در ارزیابی سیل‌خیزی حوضه است. با توجه به بررسی‌های انجام شده تجمع جریان 
از ارتفاع بالا به سمت ارتفاع Gel‏ در حوضه زیاد می‌شود و با پیش‌تر شدن تجمع جریان بر سیل‌خیزی حوضه 
افزوده می‌شود (پاتریکاکی و همکاران» ۲۰۱۸). ازاین‌جهت برای فازی‌سازی نقشه تراکم جریان از تابع خطی — 
افزاینده استفاده شد. همچنین نقشه تجمع جریان با استفاده از نقشه مدل رقومی ارتفاعی و با اکستنشن Hydrology‏ در 
نرم‌افزار Arc Gisl0.5‏ تهیه شد. 
شاخص رطوبت توپوگرافی (TWD‏ 
شاخص رطوبت توپوگرافی یک مفهوم فضایی است که انباشت جریان و توانایی حرکت آب به دلیل شیب و 
نیروی گرانش از هر حوضه را مدل می‌کند؛ بنابراین مقادیر بالا برای شاخص رطوبت توپوگرافی نشان از حضور 
مناطق مستعد با توپوگرافی مرطوب و پتانسیل تولید رواناب سطحی در حوضه است (هندینی ' و همکاران ۲۰۲۱). 
در این پژوهش برای فازی‌سازی نقشه رطوبت BUS y i‏ از تابع خطی — افزاینده استفاده شد. همچنین نقشه 
شاحص رطوبت توپوگرافی با استفاده از نقشه مدل رقومی ارتفاعی و از رابطه(۱) در Arc Gis10.5‏ به‌دست آمد. 
TWI -In(As/ Tanp) (1)‏ 
شاخص توان آبراهه (SPD‏ 
شاخحص توان آبراهه به توصیف پتانسیل فرسایشی و میزان رواناب سطحی می‌پردازده به‌طوری‌که مقادیر بالا برای 
SPI‏ بیانگر بالا بودن ظرفیت رولناب سطحی در حوضه است SL)‏ و همکاران. ۲۰۲۰). ازاین‌جهت برای 
فازی‌سازی «Y‏ توان جریان از تابع خطی- افزاینده استفاده شد. شاخحص توان جریان با استفاده از نقشه مدل رقومی 
ارتفاعی حوضه از رابطه(۲) در ag Arc Gis10.5‏ شد. 
SPI -(ASxTanp) (Y)‏ 


در رابطه‌های ۱ و AS Y‏ مساحت حوضه (m?/m)‏ و 8 شیب بر حسب درجه است (نمپاک" و همکاران ۲۰۱۶). 


انحنای سطح و انحنای طولی 


را نمی‌توان نادیده گرفت (داس " ۲۰۱۸). در انحنای سطح, مقادیر منفی» صفر و مثبت به ترتیب بیانگر مقعر» سطح 


1 Handini 
2 Nampak 
3 Das 


1۵۰ جغرافیا و مخاطرات محیطی be‏ چهارم 


زیاد و کم شدن رواناب می‌شود (کوستاج! و همکاران» ۲۰۲۱). برای فازی‌سازی لایه‌های انحنای سطح و طولی از 
تابع de‏ — کاهنده استفاده شد. انحنای سطح E‏ طولی با استفاده از آفزونه E DEM Surface‏ با استفاده از نقشه مدل 
رقومی ارتفاعی در نرم‌افز ار Arc Gis10.5‏ حاصل شد. همچنین در (Y)Error! Reference source not found.‏ 


لایه‌های_فازی شده معیارهای موثر بر وقوع سیل‌خیزی حوضه OLS‏ داده شده است. 
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شکا ۳- نفشه‌های فازی‌سازی شده معیارها 


1 Costache 


سال دوازدهم اولویت‌بندی و تحلیل فضایی پتانسیل سیل‌خیزی ... ۱۵۱ 


۳- نتایج و بحث 

۱-۳- اولویت‌بندی معیارها 

نتایج حاصل از تحلیل سلسله مراتبی و مقایسه‌های زوجی بین معیارهای انتخابی مؤثر بر سیل‌خیزی حوضه بر 
اساس نظر کارشناسان آبخیزداری OLE‏ داد که معیارهای باران» کاربری اراضی, بافت خاک» سنگ‌شناسی, فاصله از 
رودخانه شیب رتبه آبراهه تراکم زهکشیء شاخص رطوبت توپوگرافی» تجمع جریان شاحص توان جریان, انحنای 
سطح» انحنای طولی و ارتفاع به ترتیب با وزن‌های نهایی ۰/۲۲۹ ۰/۱۵۹ ۰/۱۲۰ ۰/۱۱۳ ۰/۰۷۹ ۰/۰۳۱ ۰/۱۵ 
۹ ۷ ۰۰/۰۲۵ ۰/۰۱۹ ۰/۰۱۸ و ۰/۰۱۳ بیش‌ترین تا کم‌ترین میزان تأثیرگذاری در سیل‌خیزی 
حوضه دارند شکل(؟). در پژوهش حاضر معیار باران به دلیل شروع کننده رواناب و مهم‌ترین عاملی است که به‌طور 
مستقیم در چرخحه هیدرولوژی دخللت دارد.با وزن ule‏ (۰/۲۲۹) بیش‌ترین میزان تأثیرگذاری سسیل‌خیزی را در 
حوضه دارد. در پژوهش میرموسوی و اسمعیلی (۱1۰۰) gd‏ معیار باران بیش‌ترین وزن (۰/۲۱۶) را در سیل‌خیزی 
حوضه داشت. همچنین نرخ شاخص ناسازگاری برای مقایسه‌های زوجی بین معیارها ۰/۰۸ به دست آمد که حاکی از 
قابل‌قبول بودن قضاوت‌ها است. به‌طورکلی اگر شاخص ناسازگاری کم‌تر از ۰/۱ باشد قضاوت‌های انجام شده 


سازگار و موردقبول هستند (ساعتی ؛ ۱۹۷۷). 


Priorities with respect to: 
Goal: Flood Potential 


Rain 

Land use 

Soil Texture 
Lithology 
Distance to river 
Slope 


stream order 
dragian Density 
TWI 
Flow Accumulation 
SPI 
plan curvature 
profile curvatura 
elevation 
Inconsistency — 0.08 
with 0 missing judgments. 


شکل -٤‏ مقایسه زوجی بین معیارها 


(۰/۲۲۵) بیش‌ترین تأثیر را در سیل‌خیزی حوضه دارند به‌بیان‌دیگر» کاربری‌های پهنه‌های آبی و اراضی لخت پتانسیل 


1 07 


YOY‏ جغرافیا و مخاطرات محیطی be‏ چهارم 


سیل‌خیزی بالاتری ایجاد می‌کنند (سکل۵). نتیجه پژوهش حاضر از نظر تأثیر کاربری‌های پهنه‌های آبی و اراضی 
لخت در سیل‌خیزی حوضه با پژوهش شهابی (VEO)‏ مطابقت دارد. همچنین واحدهای سنگ‌شناسی plvav‏ 
(۰/۲۵۲) و Etliv‏ (۰/۲۱۳) بیش‌ترین وزن نهایی را به خود اختصاص دادند که در سیل‌خیزی حوضه نقش موثرتری 
دارند (شسکل1). همچنین نتایج مقایسه‌های زوجی برای لایه بافت خاک از نظر کارشناسان حاکی از آن است که 
کلاس بافت خاک سنگین- خیلی سنگین و سنگین با وزن نهایی به ترتیب ۰/4۹۹ و ۰/۲۶۵ بیشترین تأثیر را در 
سیل‌خیزی حوضه دارند (شکل۷) که این امر به ماهیت نفوذپذیری این خاک‌ها در ارتباط با نفوذ آب به داخل خاک 
مرتبط است که در خاک‌های سنگین نفوذپذیری آب کم است. 
Priorities with respect to:‏ 


Goal: Flood Potential 
>Land use 


Water 
Bare Ground 
Built Area 


Flooded Vegetation 
Crops 
Grass 
Scrub/Shrub 
Trees 
Inconsistency = 0.06 
with 0 missing judgments. 


شکل 0- مقایسه زوجی بین طبقات کاربری اراضی 


Priorities with respect to: 
Goal: Flood Potential 


»Lithology 


PLvav 

Etliv 

OMviv 

URig 

Oiiv 

Qr2 

OMsm 

OMI 

Kiam! 

Inconsistency = 0.08 

with 0 missing judgments. 


شکل V‏ مقایسه زوجی بین طبقات واحدهای سنگ‌شناسی 


سال دوازدهم اولویت‌بندی و تحلیل فضایی پتانسیل سیل‌خیزی ... ۱۵۲ 


Priorities with respect to: 
Goal: Flood Potential 
»Soil Texture 


HEAVY-VERY HEAVY 
HEAVY 
MODERATE-HEAVY 
MODERATE 
MODERARE-LIGHT 
LIGHT-MODERATE 
Inconsistency = 0.08 

with 0 missing judgments. 


شکل ۷- مقایسه زوجی بین طبقات بافت خاک 


Y-Y‏ - تحلیل فضایی سیل‌خیزی حوضه 

پژوهش حاضر با هدف اولویت‌بندی و تحلیل فضایی سیل‌خیزی حوضه در مناطق نیمه‌حشک انجام شد. 
desc‏ ان ر فی ر ف با ااه ارغ مسا اشفا مار مرس عون E‏ ا bes‏ 
اس تفاده از روش FUZZY-AHP‏ در da‏ نرم‌افزار M Arc Gis10.5‏ — آمد. ces‏ حاصل از نقشه نھایی 
سیل‌خیزی به‌دست adal‏ از حوضه نشان داد که نقشه نهایی سیل‌خیزی حوضه» ارزشی بین ۰/۱۶۱-۲ دارد که به 
ون ee‏ ل کی ترس gl‏ خی E‏ ریب AN CNW =) oec‏ 
EYN) ۰/۳۶۱ -۰/۳۸۳( ۱‏ = ۰/۳۸۳) و (۱/1۳۱-۰/۵۵۳) تقسیم شد (شکل ^( همچنین در پژوهش 
حسین‌زاده 3 همکاران (AY44)‏ ارزش نهایی نفشه خحطر c»‏ در محدوده بین ۵-- ۰/۰۳۱۲ گزارش شده است که 


£ = High : 0.552 
" Low : 0.141 


M ES تا‎ S p- ~ km 
0 275 5 " 16.5 22 


5 0 230۳ 


شکل ۸- ارزش نهایی نقشه سیل‌گیری حوضه 


\OY‏ جغرافیا و مخاطرات محیطی oles‏ چهارم 


نتایج حاصل از نقشه نهایی سیل‌خیزی حوضه مطالعاتی (شکل٩) OLS‏ داد که بیش‌ترین قسمت‌های سیل‌خیز 
حوضه منطبق بر قسمت‌های مرکزی تا جنوب شرقی حوضه است که این موضوع می‌تواند به دلیل بافت خاک 
متوسط تا خیلی سنگین, وجود واحدهای سنگی Omviv Oliv Etilv‏ (به دلیل کاهش نفوذپذیری و تجمع جریان»؛ 
obe‏ باران دریافتی Y‏ ساله بیش‌تر» شیب و ارتفاع کم و وجود آبراهه با رتبه بالا باشد. در پژوهش عرب‌عامری 
و همکاران (ITA)‏ نیز بیشتر مناطق حساس به سیل‌خیزی در بخش‌های شمالی و شمال غربی واقع شده‌اند و مناطق 
مرکزی حساسیت کمی نسبت به سیل‌خیزی حوضه دارد. نتایج حاصل از مساحت‌های سیل‌خیزی (جدول۵) حاکی 
از این است که از کل مساحت حوضه حدود ۰۲۱/۰۹ YVON ۷۵/۸۷ OVAN‏ و ۳۰/۵۲ کیلومترمربع به ترتیب در 
طبقه سیل‌خیزی خیلی کم. کم. متوسط زیاد و خیلی زیاد قرار دارند. به‌این ترتیب بیش‌تر مساحت حوضه در طبقه با 
سیل‌خیزی متوسط تا خیلی زیاد قرار دارند؛ بنابراین Lb‏ اقدامات حفاظتی و JES‏ سیل در این حوضه موردتوجه 
قرار گیرد. همچنین ۲۷/۰۵۳ ۱۲/۶۰7 درصد حوضه به ترتیب از پتانسیل سیل‌خیزی زیاد تا خیلی زیاد برخوردار 
است. در پژوهش داس (۱۰۱۸) نیز ۲۰ درصد حوضه در طبقه بسیار زیاد گزارش شده است. از دیگر نتایج این 
پژوهش می‌توان به این موضوع اشاره کرد که مناطق با سیل‌خیزی خیلی کم تا کم منطق بر قسمت‌های انتهایی حوضه 
است که می‌تواند ناشی از عوامل همچون فاصله زیاد تا آبراهه وجود سازندهای زمین‌شناسی 012 (نهشته‌های 
تراسی کواترنر» میزان بارش ۳۱ ساله دریافتی کم‌تر و جنس بافت خاک سبک تا متوسط a)‏ دلیل نفوذپذیری بالا) و 
تمرکز بالای آبراهه‌های با رتبه ۱ اشاره کرد. همچنین روستاهای قمصر و قزاآن در پهنه سیل‌گیری زیاد تا خیلی زیاد 
قرار دارند لذا اقدامات مدیریتی - حفاظتی برای این روستاها و اراضی کشاورزی و باغ‌های آن‌ها امری ضروری است. 


20 Kilometers 


rece 91 23001 


شکل -٩‏ نقشه نهایی تحلیل فضایی سیل گیری حوضه 


سال دوازدهم اولویت‌بندی و تحلیل فضایی پتانسیل سیل‌خیزی ... 100 


جدول ۵- مساحت و درصد سیل‌گیری حوضه در هر طبقه 


طبقه‌بندی سیل‌گیری مساحت (km?)‏ درصد 
خیلی کم ۳۱/۳۹ NOVA‏ 

کم 041 YAYI‏ 
متوسط Vo/AV‏ ۱۳۰/۸۶۱ 
35 1701 7/۰07 
خیلی زیاد ۳/۹۲ WES‏ 


-٤‏ نتیجه گیری 

سیلاب یکی از مخاطرات محیطی است که با ایجاد خسارت‌های اقتصادی» اجتماعی و زیست‌محیطی منجر به از 
54 0555 6.05 شنده که lee cele a‏ ان از ساین acd eles‏ کی ر yl poa‏ اا اوا ux‏ 
سابقه سیل در منطقه مورد‌طالعه طی سال‌های اخیر تحلیل فضایی سیل‌خیزی حوضه و شناسایی مناطق با سیل‌خیزی 
بالا از اهمیت ویژه‌ای برخوردار است. این پژوهش با اهداف شناسایی و اولویت‌بندی معیارهای phe‏ بر سیل‌خیزی 
حوضه و تحلیل فضایی سیل‌خیزی حوضه انجام شد. بدین منظور در راستای هدف پژوهش حاض از روش 
FUZZY- AHP‏ برای اولویت‌بندی و تحلیل پتانسیل سیل‌خیزی حوزه آبخیز قمصر از VE‏ معیار (لایه‌ها) «QUU‏ 
ارتفاع کاربری اراضیی بافت خاک. سنگ‌شناسی, فاصله از رودخانه شیب تراکم زهکشی رتبه آبراهه» شاخص 
رطوبت توپوگرافی. شاخص توان آبراهه» تجمع جریان» انحنای سطح و انحنای طولی استفاده شد. نتایج حاصل از 
تحلیل سلسله مراتبی و مقایسه‌های زوجی نشان داد که معیار باران با وزن نهایی ۰/۲۲۹ مهم‌ترین عامل موثر در تحلیل 
پتانسیل سیل‌خیزی حوضه است. نتایج حاصل از نقشه نهایی سیل‌خیزی حوضه نشان می‌دهد. از کل مساحت 
حوضه حدود ۷۵/۸۷ کیلومترمربع از پتانیل سیل‌خیزی متوسط برخوردار است که با توجه به Sep)‏ این مناطق از 
قابلیت تبدیلی بالایی به مناطق با سیل‌خیزی زیاد تا خیلی زیاد برخوردارند؛ لذا این مناطق باید برای اقدامات مدیریتی 
کنترل سیلاب و آبخیزداری در اولویت قرار گیرند. همچنین بخش‌های مرکزی و امتداد جنوب شرقی حوضه در طبقه 
سیل‌خیزی زیاد تا خیلی زیاد قرار دارند که این می‌تواند ناشی از میزان باران دریافتی بیش‌تر, بافت خاک متوسط تا 
خیلی سنگین, سنگ‌شناسی نفوذناپذیر Omviv Oliv Btilv‏ شیب و ارتفاع نسبتاً کم تمرکز آبراهه با رتبه بالا (۳ و 
(E‏ در این بخش‌ها باشد. همچنین دو روستای قمصر و قزاآن در محدوده سیل‌خیزی زياد تا خیلی زياد قرار دارند که 
سیل می‌تواند به اراضی بافت روستایی و کشاورزی و باغ‌های این روستاها خسارت وارد کند. ازاین‌جهت AL‏ 
اقدامات حفاظتی برای این روستاها موردتوجه قرار گیرد. با توجه به نتایج پژوهش حاضر می‌توان نتیجه گرفت که 


10۶ جغرافیا و مخاطرات محیطی be‏ چهارم 


سیل‌خیزی بالا حوضه و ارائه برنامه‌ریزی‌های مدیریتی برای منطقه به همراه داشته باشد؛ cyl ple‏ آگاهی از پتانسیل 
سیل‌خیزی حوضه می‌تواند برای تدوین برنامه‌های Ce pte‏ بحران در هنگام مواجه‌شدن با سیلاب مؤثر واقع شود. با 
bas‏ اثرات آن در محدوده دارای طبقه‌بندی سیل‌خیزی متوسط تا خیلی ob‏ از ساخت‌وساز اصولی استفاده گردد و 


سازه‌هایی که حساس به تخریب در اثر سیلاب هستند مقاوم‌سازی صورت ae‏ 
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